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SUMMARY 


This is a working paper in which a formulation is given for solving the boundary-layer equations in 
general body-fitted curvilinear coordinates while retaining the original Cartesian dependent variables. The 
solution procedure does not require that any of the coordinates be orthogonal, and much of the software 
developed for many Navier-Stokes schemes can be readily used. A limited number of calculations have 
been undertaken to validate the approach. 


INTRODUCTION 


The boundary-layer approximation is a useful engineering tool which has contributed significantly 
to the understanding of viscous flow at high Reynolds number. The boundary-layer equations require the 
use of a body conforming coordinate system and the flow Reynolds number must be high. In developing 
the usual boundary-layer equations, both the independent variables and the dependent velocity variables 
are transformed to the new body conforming coordinates. For body surfaces with little curvature, the 
boundary-layer equations cast in terms of the new dependent variables more or less simplify back to a 
flat plate or Cartesian-like-form of the equations along a developed surface. However, if the body has 
appreciable curvature, the equations become more complicated. They are particularly more complex if a 
nonorthogonal coordinate system is used, yet for many applications it is difficult to generate an orthogonal 
coordinate system along the body surface. 

In this note, a formulation for the boundary-layer equations in terms of the original Cartesian velocity 
variables is described for body-fitted general curvilinear coordinates. Although the collaborating compu- 
tational experiments that have been undertaken are limited, the proposed alternate form of the governing 
equations may offer several advantages in terms of numerical stability by avoiding coordinate source terms. 
Moreover, this alternate form of the boundary-layer equations does not require that any of the coordinates 
be orthogonal, and software (grids, boundary condition routines, etc.) developed for many Navier-Stokes 
schemes can be readily used. 

This formulation was partially motivated by discussions with Dr. H. Yoshihara. 


BACKGROUND 


The compressible boundary-layer equations for the unsteady, three-dimensional flow of a perfect gas 
over a flat plate can be written 


Pt + ( pu)x + ( pv)y + ( pw) z = 0 


(la) 


pU t + puu x + pVUy + pWU z - -p x + ( pu y )y 


( 16 ) 


1 



( 1 C) 


PWt + pUWx + pVWy + pWW z = —p z + (p,Wy)y 


pn t + pun a 




. Pr L 1 


Pr — 1 2 2 s 

+ ~ (u +w 2 ) j 


■Jh 


+ p 


(Id) 


where p and p are the density and pressure, u , v , w are Cartesian velocity components, and H is the specific 
total enthalpy, H = “j 2 with e the total energy per unit volume. Here Pr is the effective Prandtl number. 
The equations are nondimensionalized and y is a coordinate normal to the body surface. These equations 
can be used for bodies of slight surface curvature using x and z as distances along the surface with u, w, 
and v the corresponding velocities. 

If surface curvature effects are taken into account, the boundary-layer equations take on a more com- 
plex forni primarily because of the addition of coordinate source terms (refs. 1 and 2; Warsi, Z. U. A.; 
unpublished notes, 1984). For example, if f , rj, and £ are defined as shown in figure 1, with p normal to 
the surface, Panaras (ref. 2) derived the form 

(/-'p)r + ( J~'pU)( + U~ X pV\ + (J~ x pW) ( = 0 (2 a) 


pU t + pUU( + pVU v + pWU ( + K ai pU 2 + K a2 pUW + K a3 pW 2 = 

^a5 Pf ^ + PiK a $ JVji 

+ Jd„[-^j(U v + K a6 u + K a7 W)-i (26) 


pWt + pUW( + pVW v + pWW$ + Kb\pU 2 + KbipUW + KbipW 2 = 

Kb4P( — Kbspc + pKb% U v + pKb9 

+ Jd v [-£y(W v + Kk,U+ K bl W)] (2c) 

922 J 


( /- 1 e) T + [ J- 1 ( e + p) !7] f + [ /"> ( c + p) V], + [ J~ l ( e + p) W\ = 


d v [pK e id v U 2 + pK e2 d n (UW) + pK €i dJY 2 + 


P 

Pr(7 - 1)022.1 



(2d) 


where the metrics are given in appendix 1 and U and W are contravaiiant velocities as defined later. Here 
the 77 -coordinate must be normal to the body surface. 

Although there are more terms than previously, this is a fairly clean form of the boundary-layer equa- 
tions. From a numerical point of view, however, equation (2) are more difficult to solve than equation (1) 
because derivatives of metrics must be formed which are not always smooth and because the extra source 
terms can adversely effect numerical stability. For example, any coordinate source term such as K a \pU 2 
in equation (2b) can degrade stability, in this case if K a 1 U has a negative value. 
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In the following section the steps to bring the Navier-Stokes equations in Cartesian coordinates into 
general curvilinear coordinates and velocity components are begun as if to derive equations like equation 
(2). Here these steps are only taken to determine how to group the equations. Once the pressure in the 
viscous layer is determined from the uncoupled normal-like momentum equation, primitive forms of the 
momentum equations are used in place of equations 2b and 2c. The primitive forms with given pressure 
and use of a thin-layer approximation are easier to treat numerically than equation ( 2 ). 


DEVELOPMENT OF CURVILINEAR BOUNDARY-LAYER EQUATIONS 


The three-dimensional Navier-Stokes equations in Cartesian coordinates can be written as 


(p)t + (pw)i + (pv) y + (pw) z - 0 (3a) 

pu t + puu x + pvu y + pwu z + p x = R xmom ( 3 6) 

pv t + puv x + pvvy + pwv z + p y = R ymom ( 3 c) 

pw t + puw x + pVWy + pww z + p z = Rzmom ( 3 d) 

pHt + puH x + pvHy + pwH z -Pt = Revet ( 3 e) 

where ( Rxmom > Rymom ) Rzmom ) Revet will represent the viscous terms. 

Transforming the independent variables x, y, and z to body conforming coordinates ( , rj, and ( gives 

( J -1 p) T + ( / pU\ + ( J~ l pV\ + ( J~ l pW\ = 0 ( 4 o) 

pu t + pUu( + pV u v + pWu$ + (6^ + r) x p v + C xpO = R xmom (46) 

P v t + pUV£ + pVUr, + pWv C + (tyPz + TJyPfj + £ yP< ) = Ry mom (4 c) 

pw t + pUw i + pVw v + pWw C + (&p e + n z p v + Cpc) = Rzmom (4 d) 

pHt + puHl + pvH v + pwH{ - Pt = Revet ( 4 e) 

where the U, V, and W represent unsealed contravariant velocities; e.g., 

U - £ x u + (yV + £ z w (5 a) 

v = VxU + 7 ] y v + r\ z w ( 56 ) 

W = CxU + (yV + ( 2 w (5c) 

and J is the transform Jacobian 

J = (xty v z<; + x^z v + x v y<;Z( - x^z v - x v y£Z$ - a^y ^)” 1 


and other metrics are defined at the end of appendix 1. Note that the Cartesian velocity variables u, v, 
and w are still retained as dependent variables and that the momentum equations are still the Cartesian 
momentum equations. 
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To obtain equations such as equation 2 using new dependent velocity variables U, V , and W, the 
dependent velocity variables can be transformed by taking linear combinations of the Cartesian momentum 
equations. Specifically these equations can be multiplied by the matrix C defined as 

€x Cy 

C = §* fjy f} 2 (6) 

- C* C y C z - 

where the overbar implies scaling of the metrics such that ( x = ix/^Jil + fj x = 

W + Vy + Vz etc. Multiplying the three momentum equations by C gives 

pUt + pUU( + pVU v + pWU $ + + Vf • (Vfp^ + Vtjpj, + VC P() = V( ■ R (la) 


pV t + pUV£ + pVV v + pWV$ + s v + Vtj • ( V£p£ + Vrjp v + VCpc) = Vp • R (lb) 

pWt + pUW{ + pVW v + pWW c + s < + VC * (' VCp^ + Vpp n + VCp<) = VC • R (7c) 

where R = ( R xm om, Rymom , RzmomY and where V is the gradient operator. Terms like Vf • VCpc should 

be interpreted as (f«C* + C»(» + f.C*)Pc- 

The terms s^, s v , and are coordinate source terms given as 

$Z = - pU(ud(l x + vd^y + wd(C z ) 

- pV ( ud v £ x + vd v (y + wd v ( z ) ( 8 o) 

- pW ( ud^ x + vdtfy + wd{£ 2 ) 

Srj = — pU(ud{fi x + vd(fj y + wd{fiz) 

- pV ( ud v fi x + vd v fjy + wd v fjz) ( 8 b) 

- pW ( ud^fjx + vd^fjy + wdvfjz) 

•s< = - pU(ud( Cx + vd{( y + wd(Cz) 

- pV ( ud v ( x + vd v (y + wd v C) ( 8 c) 

~ pW ( ud^Cx + vd£ y + wdc£z) 

The Cartesian velocities in these source terms can be replaced with contravariant velocities using 

( u \ [a* x v rr c 1 /U\ 

[ v I = i/e ft, J/c V (9) 

\w/ Lz£ z v z^WwJ 

where the transformation matrix is simply O” 1 for C unsealed. Often the source terms are rearranged and 
written with Christoffel symbols, but such notation is not needed here. 
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In the body conforming coordinates, the wall corresponds to 77 = const, so V77 is normal to the wall 
and the flow tangency condition requires that V = V17 • q - 0 where q = ( it, v, w)*. For a thin viscous 
layer near the wall, V will be 0 ( 6 ) (as will the variation of ( and ( metrics with 77, i.e. d v £ x = 0 ( 6 ) etc.). 
For small V the 77-momentum, equation 7 b, simplifies to 


5tj + V77 • ( Vfp£ + Vrjptj + VCpc) = Vtj • R 

Since R represents viscous stress terms along the wall, the product V77 • R should be zero at the wall and 
small away from the wall. Consequently the 77-momentum equation furthers simplifies to 

s n + V77 ■ ( vfpf + Vr)p v + VCpc) - 0 ( io) 

Given outer-edge conditions for p, this represents a simple wave equation with a source that can be used 
to evaluate the pressure through the viscous layer. This equation further simplifies if the 77 coordinate is 
orthogonal to the surface as V77 • Vf = 0 and Vtj • VC = 0 . Even for nonorthogonal coordinates, the 
coefficient V77 • V77 is generally so much larger than the other coefficients that the assumption 

Pv = 0 (H) 

is often valid over a thin layer near the wall. Thus from equation ( 11 ) p is prescribed throughout the 
boundary layer to its specified edge value. Therefore, pressure does not depend on the other dependent 
variables within the boundary layer if equation (11) is used, while the dependency is weak if equation (10) 
is used. 

Once the pressure is determined in the boundary layer from the 77 momentum equation, the pres- 
sure derivatives in the Cartesian momentum equations, equations 4 b-d can be evaluated. With pressure 
specified, the inviscid portion of the Cartesian momentum equations using only transformed independent 
variables are very easy to solve for u, v, and w. The inviscid equations with specified pressure are just 
simple convection equations. Consequently, assuming the thin-layer viscous terms cause no difficulties, 
equations 4 b-d can be readily solved for it, v , and w in place of the more complex equations 7 a and 7 c for 
U and W. Taking a linear combination of (u, v, wY with Vf and VC (now ignoring metric scaling), the 
transformed velocities U and W are then readily formed as 

W = VC • q 

At this point the same linear combination of Cartesian momentum equations is used to predict U and W , 
but the Cartesian momentum equations rather than the ( and C momentum equations are used directly. As 
a consequence the complex source terms are avoided. 

Up to this point the viscous terms have been mostly ignored in the development. Their complexity 
remains to be checked, as well as whether they tend to couple the equations together. 

To examine the viscous terms, it is illustrative to first consider incompressible flow. For incompress- 
ible flow in which the coefficient of viscosity is constant the viscous terms are given as 


2 

Rxmom ~ /lV It 

(12 a) 

Rymom = M V V 

( 126 ) 

Rzmom ~ /lV W 

(12c) 
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Clearly these terms do not couple the momentum equations together. Now the Laplacian operator trans- 
forms (in divergence form) as 

v 2 = v$ • vod( + dtj-'ivt ■ Vv)d v + a { j-'(v$ • vofy 

+ V'(Vi j • VOd( + d„J-'(Vv ■ Vt))3„ + d v J-'(Vv ■ VQd ( (13) 

+ afcj-'cvc • voa f + vHvc • v V )d v + d ( j~\v c • voa<] 

The extra cross derivative terms caused by coordinate transformation have been encountered in potential 
and Navier-Stokes codes and can be differenced so as not to cause numerical instabilities. However, equa- 
tion (13) can be further simplified by making a thin-layer approximation in which case these terms can be 
eliminated altogether. Dropping all derivatives with respect to ( and the Laplacian is reduced to 


V 2 = Vtj • Vrj)9,,] (14) 

which is particularly simple. Thus for an incompressible viscous term and specified pressure, the momen- 
tum equations (4b-d) are uncoupled and are readily solved for u, v, and w. 

If the Cartesian form of the compressible viscous terms undergoes independent variable transforma- 
tion, ( x, y, z) to (£, 77, 0 , and if they are subject to a thin-layer approximation, they simplify to 


Rxmom = JdriiJ '(limiUr) + (/i/3) 9722 Vx)] (15o) 

Rymom (/^rUiU^ + (/^/3) Tfl2 tjy) ] ( 156) 

Rzmom = Jd v [ J~ l (lim\W v + (^/3)77t 2 77*)] ( 15c) 

Rener = Jd v {J~ l [ fJ,m\ m3 + (/i/3)m 2 ( T} X U + TfyV + r} z w)]} (15 d) 


where mi = vl + vl rn 2 = and m 3 = [(u 2 + v 2 +w 2 ) v ]/2+Pr~ l (l-l)~ l (a 2 )r ) . 

and J is again the transformation Jacobian, Here the viscous terms were first put into divergence form, and 
then simplified. As a result they are identical to the viscous terms used in many thin-layer Navier-Stokes 
codes. This can cause some error for bodies with high curvature as metric terms are also being discarded. 
Making the thin-layer approximation on the nonconservative form of the viscous terms eliminates this 
problem. This form is given in appendix 2. 

The compressible viscous terms are coupled through the m 2 term, but the coupling appears to be 
weak and has been treated explicitly in many numerical schemes without degradation of the time step or 
iterative convergence. Moreover, because the coupling occurs in only one direction, rj, it would not be too 
costly to account for it in a fully implicit manner. 


CURVILINEAR BOUNDARY-LAYER EQUATIONS 


In summary, a form of the boundary-layer equations for general curvilinear coordinates is given by 
(with rj chosen away from the surface): 
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Normal Momentum 


Pv = 0 


£ and C Momentum 

pu t + pUu ( + pV u n + pW U< + (£ xP( + Vx p v + CxP<) = Jd v {J- x lpmiu v + (p/3 )m 2 p x ]} 


pv t + pUv ( + pVv n + pWv( + (( yP( + T) y p v + ( yP< ) = Jd^J-'lpmiv,, + (p/3 )m 2 n v ]} 

pw t + pUw i + pVw„ + pWw ( + (&p £ + t) 2 p„ + ( zP< ) = Jd v {J-'[p,mvw„ + (p/3 )m 2 tj 2 ]} 
where two linear combinations of u, v, and w are used 

U = Vt-q 

W = V£ ■ q 


Energy 


pH t + puH{ + pvH v + pwH ^ — pt - JB^J 1 [pm\m ,3 + (/x/3) m 2 (rj x u + P y v + rj 2 w)]} 
where 

™i = nl + nj + v, 

m2 = rjxUf, + TjyVq + 7] z Wfj 

m 3 = [(u 2 + v 2 + tn 2 ),]/2+ Fr _1 ( 7 - ^-'(o 2 ), 

Constitutive 


P _ pToo 
poo 'Epoo 

where 


r_(T-D rrr (u 2 + v 2 + w 2 ). 
Toe al 2 J 

Continuity 


U 'p)r+ (J-'pU)( + (J~ l P V\ + (J~ l pW)( = 0 


( 11 ) 

(46) 

(4c) 

(4d) 

(4e) 

( 16 o) 
(166) 
(4 a) 
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The six equations, normal momentum (eq. 11), two linear combinations of (eq. 4b), (4c), and (4d) 
to form U and W , energy (eq. 4e), state (eq. 16a), and continuity (4a), can be used to determine the six 
variables p,U,W, H, p, and V. The Cartesian velocity components are then obtained from 



Generally p is determined from 





*<“ 


V( 

Vv 

y< 



z v 

z O 


Pri 

= 0 




(9) 


( 11 ) 


and an outer-edge boundary condition so that p = p e dg e along r\ coordinates. In more general cases pressure 
can be determined from 

fir, + Vrj • ( V£p* + Vr)p v + VCp<) - 0 ( 10 ) 

using the same outer-edge condition with the source term defined as 


St, = — pU(ud(f} x + vd(f} y + wd^) 

- pV ( ud n rj x + vd n fj y + wd v f ) z ) ( 8 b ) 

— pW(ud(f) x + vd(fj p + wd^riz) 

In this case pressure will vary throughout the boundary layer. 

For C (instead of rj) away from the surface, linear combinations of U and V are used, p c = 0 , mi = 

Cx + Cy + Cz » m 2 = C* u < + CyV( + ( Z W(, m 3 = [(u 2 + v 2 + w 2 \] / 2+ Pr -1 (7 — l) _ 1 (a 2 )^, and viscous 
derivatives are taken with respect to 


NUMERICAL TESTING 


This formulation has been tested in steady state applications using the time-like boundary-layer 
scheme reported by Van Dalsem and Steger in (ref. 3). For a prescribed edge pressure, the equations 
are solved in the following way with p« = 0 . Using central differencing in rj and upwind differencing in 
£ and equations (4b) to (4e) are used to update u, v, and w, and H. As pressure was already obtained 
using one linear combination of the momentum equations, only two linear combinations of the momentum 
equations can be used to determine the velocities. Thus, U and W are formed from the momentum updated 
w ' The equation of state, (16a), is used to update p, with T defined from equation (16b). The 
third linear combination of u, v, and w is obtained by solving continuity for V using already updated U, 

W, ana p The final updated form of the Cartesian velocity components are then obtained from updated 
U , V , and W using equation ( 9 ). 


Two test cases were used to verify the algorithm. The first case, flow over a flat plate, was simply 
used to verify that the algorithm was coded correctly. To bring in three-dimensional effects and to verify 
some or the metric terms, the computational grid was rotated on the flat plate as sketched in figure 2. In this 
simple test, a compressible Blasius solution was specified at all boundaries. Figure 3 shows the Blasius 
result and a typical computed profile (with Pr = 1) taken from the center of the 20 x 10 x 30 grid for a 
grid that was not rotated. Significant relative error is detected for the vertical velocity, but the overall error 
^small. Figure 4 shows Cartesian velocity profiles for a similar calculation but with the grid rotated by 
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A more interesting test case was provided by using the boundary-layer algorithm to verify a computed 
thin-layer Navier-Stokes result on a prolate spheroid (ref. 4). In this case the Navier-Stokes result was first 
obtained for the body at 10 ° angle of attack, Moo -0.17, and Re - 7 .7 x 10 6 based on the body length. 
In the experiment (ref. 5) the boundary layer was tripped at the x/L ~ 20% station. In the calculation 
the trip was somewhat taken into account by using laminar values of the viscosity coefficient over the first 
20% portion of the body and using a eddy viscosity turbulence model from that point on. Having obtained 
the Navier-Stokes result, the new boundary-layer code was also run on a portion of the Navier-Stokes grid 
from x/L = 5% to x/L = 80% . Starting boundary-layer profiles and edge conditions were taken from the 
Navier-Stokes calculation, the “edge” corresponding to the 25 th grid point up from the wall of the Navier- 
Stokes grid (about 0 .003 body lengths). Various profiles from this calculation are shown in figures 5 and 6 
at x/L = 0 .48 and x/l = 0 .64 . Also shown in figure 7 are boundary-layer-computed wall-turning angles, 
Navier-Stokes-computed wall-turning angles, and experimental data (ref. 5). Although discrepancies are 
evident, the overall agreement is very good on the windward portion of the body and is more than adequate 
to verify the boundary-layer equations and algorithm. 


CONCLUDING REMARKS 


A formulation has been given for the boundary-layer equations using general body-fitted curvilinear 
coordinates and retaining the original Cartesian dependent variables so that coordinate source terms are 
avoided. The formulation does not require that any of the coordinates be orthogonal and gridding and 
software developed for many Navier-Stokes schemes can be readily used. 

The curvilinear boundary-layer equations given previously have obvious similarity to the thin layer 
Navier-Stokes equations. However, in this boundary layer formulation, while three momentum equations 
are solved for three Cartesian velocity components, only two linear combinations of velocity variables are 
actually taken from momentum. A third linear combination of the momentum equations is used to provide 
the simplified 77 -momentum equation with a prescribed edge pressure. Pressure throughout the boundary 
layer thus is specified and uncoupled (or weakly coupled if eq. ( 10 ) is used over eq. ( 11 )) from the 
other dependent variables within the boundary layer. Tliis then allows the use of boundary-layer solution 
algorithms. 
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APPENDIX 1 


METRICS 


The metric relations needed for the Panaras curvilinear boundary layer equations 


K a i = 0 5g n %!- + g 13 ^. - 0.5g 13 9dn 


ae 


dt 




rr _ „n Sgn 13 9g33 

•K a 2 - 9 - 57 - + g 




dt 


K a 3=0.5g I3 ^ + s »^-0.5g" 9933 


9C ’ d( 

922913 




iT a4 = 

£a5 = 


0 

922933 


Ka6=0.5g"^ + 0.5g' 3d9u 


K a 1 = 0.50 


dr\ 
11 ^011 


dr) 

Kai = g 22 K a (, 


drj 


+ 0.5g nd9M 


drj 


K a9 =g 22 K al 


Kb\ - 0.5 g 


13 dgn 


+ 0 


33 


#013 


-0.50 


33 


dgn 

dC 


K b2 = g 13 9911 - - 33 9?33 


3C 


+ 9 " 


dt 


^3 = 0.5g 33 ^ + g’ 3 ^-0.5g 139933 


ac ' ' 

922 013 


3f 


^64 = 

Kb 5 = 


0 

011022 


are given as 
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K b 6 = 0V s ^ + 0.5g 33 ^ 
Kv, =0.5g ,3 ^2. + 0.5 9 33 ^ 

01 7 07} 

Kb% = Q n Kte 


tr „22 tv - 

A&9 = 0 itfi7 


*el = 


011 

2022 


if e2 = — 
0 22 


tf e 3 = 


033 

2022 


011 = O^ICf + + Z£Z{ 

022 = XfjXfj + 

033 = X(X( + + Z£Z{ 

012 “ %£Xij + + 2£2tj 

013 = X£X$ + + Z£Z{ 

023 = XnX{ + y v y{ + z v z < 


0 H = (022033 023023) / 0 

0 22 = (011033 - f?13013) / 0 

0 = (022011 — 012 012) /0 

0 = (023031 — 021033)/0 

0 = (021032 “ 022 013 ) /g 

23 

0 = (021013 - 911923) /g 

0 = (v7r 1 

The metrics used throughout the paper are defined as 

fx = J(y v zc-y<z v ) 
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6y = J ( x c z v ~ x v z 0 

6 = J( x tiV$ — x ^Ptj) 

Vx - J(y^ - y^) 

Vy = J(X(Z{-X{Zt) 
riz - J (x<;yt - x^) 

Cx = J (ytzr, - y^) 

C V = J(x n Z£ — X(Z V ) 

Cz = J( x (yv~ x vVO 

with 

= ( x tvv z < + x cyt z v + x vy< z t - x ty< z n ~ x nvt z c - x <vv z O 
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APPENDIX 2 


THIN-LAYER VISCOUS TERMS 


Let the Cartesian form of the compressible viscous terms undergo independent variable transforma- 
tion, (x,y,z) to (£, 77 , 0 in chain rule conservative form. If they are subject to a thin layer approximation 
with £ taken as the direction away from the body, they simplify to 


Rxmom ~~ 

Cx0clXf*U< + Crti<) + + C*W<)] 

+Cy^<tMCf + CxVc)] 

+C*9c[./i(Cittc + C*wc)] 


■ymom ~ 

Cxd^[ + Cx v <)] 

+Cy^c[M(CyV< + Ci;^) + (Cx^C + Cp v < + C*^)] 
+Czddt J, ( < > z v < + Cy^c)] 


'zmom ~ 

Gxfy[/J'(CzU'( ■*" Cx™<)] 

+Gyfy[v(G* v < + Gv w <)] 

+Gzd{Ui,(Gzw<: + GzwO + MC*u< + Oc + C^c)] 


C x d{[KPr *(7 - 1 ) ^(a 2 )^ + u\(G x U( + + Gz^O 

+ U^(Cx«C + Cxtt<) + V/i(Cy«< + GxV{) + WfiiGzUC + CxW<)] 

+Cy^ct^' Pr_1 (7 - l)- 1 Cy(a 2 )c + vA(Cx«c + Cy v < + 

+u/x(Cyii< + Cx^c) + v MCy v C + GyV() + c + Cy^c)] 

+G^ct«Pr “ 1 ( 7 - l) _ 1 C(a 2 ) c + w\(GxH + GyV( + GzW '<) 

+ U/i(C 2 U^ + G X W() + Vfl(G z V{ + (yW() + Wfl(GzW^ + £,«/()] 

where X = —( 2 / 3 )/i. Simply replace £ with 77 if 77 is to be used as the direction away from the body 
surface. 
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Figure 4.- Computed velocity profiles for u, v, and w taken at the center of a small grid compared with the 
Blasius solution; grid rotated by 30 ° and Pr = 1 . 







Figure 5.- Tranformed and physical plane computed boundary-layer velocity profiles compared to com- 
puted thin-layer Navier-Stokes profiles for a 10°-angle-of-attack prolate spheroid at station x/L = 0.48 
and circumferential locations of a) </> = 0°, b) $ - 60°. 
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Figure 5.- Concluded, c) <f> - 120°, and d) (f > » 180° (leeward). 
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Figure 7.- Wall- turning angles versus circumferential angle (f> for the boundary-layer and thin-layer Navier- 
Stokes computations compared with experiment for a 10 °-angle-of-attack prolate spheroid at axial stations 
a) x/L = 0 ,48. 
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y, SHEAR STRESS ANGLE 7 , SHEAR STRESS ANGLE 
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Figure 7.- Concluded, b) x/L = 0 .65 . 
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